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Abstract 

We explore the phenomenon of emergent Lorentz invariance in strongly coupled 
theories. The strong dynamics is handled using the gauge/gravity correspondence. We 
analyze how the renormalization group flow towards Lorentz invariance is reflected in 
the two-point functions of local operators and in the dispersion relations of the bound 
states. The deviations of these observables from the relativistic form at low energies 
are found to be power-law suppressed by the ratio of the infrared and ultraviolet scales. 
We show that in a certain subclass of models the velocities of the light bound states 
stay close to the emergent 'speed of light' even at high energies. We comment on the 
implications of our results for particle physics and condensed matter. 



1 Introduction 

Lorentz invariance (LI) is one of the best tested symmetries in nature [1]. From this fact 
one usually infers that Lorentz symmetry is simply a fundamental property underlying the 
particle physics and gravity. In this paper, we will entertain the opposite idea: is it possible 
that LI is not a fundamental but only an accidental symmetry of the low energy world ? 
Can LI be an 'emergent symmetry' ? 
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One motivation to insist on fundamental physics possibly being Lorentz-violating (LV) 
may come from the recently discovered consistent non-relativistic models of gravity [2, 3], 
which could by themselves be UV complete. But we will try not to attach too much to 
any concrete model and simply assume that at energies E > M for some high energy scale 
M physics is Lorentz violating. The question that one asks then is: is it possible that the 
deviations from LI become small al low energies ? If this were to happen naturally {i.e., 
without fine-tunings of parameters) it would represent an example of an emergent Lorentz 
invar iance. 

In fact, it has been known for a long time that this phenomenon does happen quite 
generically. Nielsen et al. considered in [4, 5] what happens to a LI quantum field theory 
when it is perturbed by LV operators. Let us denote collectively the couplings in front of 
these operators klv- One can then compute the renormalization group (RG) flow of these 
couplings and finds that they vanish towards the IR. The exercise has been repeated in a 
number of contexts and with diverse field contents [6, 7, 8] giving always the same result: 
the Klv A° w to zero and so in the IR the theory is attracted towards LI. Moreover, the same 
phenomenon has been observed in the condensed matter context where, in certain materials, 
effective LI emerges from intrinsically non-relativistic Hamiltonians [9, 10, 11]. 

These findings agree with the following general argument [12, 13], which shows that under 
very broad assumptions LI fixed points are IR attractive and thus quite generically can serve 
as the end points of the RG flows. Consider a LI scale invariant theory (SFT) in d space-time 
dimensions. It is believed that if this theory is unitary, it enjoys a larger symmetry, namely 
the full group of conformal transformations, i.e. it is a conformal field theory (CFT), see 
[14, 15]. Its deformation away from LI can be equivalently understood as a perturbation of 
the original Lagrangian by primary operators 1 with free Lorentz indices; schematically, 

C-CFT ^ £>CFT + K-LV^fiL.-fii ■ 

It is convenient to write the perturbation in the Lorentz-covariant form by contracting the 
free indices with appropriately chosen spurion fields, which can be classified according to the 
irreducible representations of the Lorentz group (cf. [16]). Let us for simplicity concentrate 
on the case of deformations preserving the symmetry under purely spatial rotations. Then 
all relevant spurions can be reduced to the single time-like vector u^, such that u° — 1, 
u % = 0, i = 1, 2, . . . , d — 1, and the deformation takes the form, 

CcFT ^ CcFT + KlvO^.^u" 1 . . . U» l . (1) 

In other words, it is given by the time-components of traceless symmetric tensors. Now, in 
unitary CFT's the dimensions of such operators are bounded from below [17, 18, 19], 

dim<V.. w > d-2 + l , 

where / is the number of Lorentz indices carried by the operator. Thus, all operators with 
more than two indices are automatically irrelevant. Among operators with two indices in 

1 Or so-called quasi-primaries in the case d = 2. Descendant operators are total derivatives and adding 
them to the Lagrangian does not have any effect. 
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general there is only one — the energy-momentum tensor — which is exactly marginal, 
dimTjj,, = d, the others being irrelevant 2 . However, the contribution of T^ u in (1) can be 
absorbed by the rescaling of the time coordinate [13]. Thus, the only dangerous perturba- 
tions, i.e. those that can drive the theory away from the LI fixed point, are related to vector 
operators with the dimensions in the window d — 1 < dim M < d. If no such operators exist 
in the theory, the LI point is stable. Alternatively, the vector deformations can be forbidden 
by a discrete symmetry, such as T-parity or CPT 3 . 

To avoid misunderstanding, let us stress that the above argument does not imply that 
any RG flow must end up at a LI fixed point. This is clearly not true for non-relativistic 
SFT's with dynamical critical exponent z ^ 1. What this argument suggests instead is 
that the emergence of LI at low energies is not something exceptional and must be rather 
common. 

This answers the question formulated at the beginning of the paper basically in the 
affirmative. However, in practice there is one important qualification. In weakly coupled 
theories, there are in general LV operators that are marginal at the tree level. Important 
representatives of this class are the operators describing the differences in the maximal 
propagation velocities of various particle species [21], which within the Standard Model are 
tightly constrained by experiment. Upon acquiring positive anomalous dimensions due to 
quantum corrections the coefficients of these operators do flow to zero in the IR, but too 
slowly. The K^y run logarithmically with the RG scale and for a fixed amount of running the 
suppression of the LV couplings in the IR is not sufficient to fulfill the experimental bounds. 

There are at least two possible ways to overcome this difficulty: 

(i) Endow the field theory with enough symmetries in order to reduce the number of 
allowed operators of lower dimensions, so that LI becomes an accidental symmetry of 
the renormalizable (i.e. composed of operators with canonical dimensions up to 4) 
part of the Lagrangian. This idea is realized in the mechanism by Groot Nibbelink 
and Pospelov [22], based on a non-relativistic form of supersymmetry. Combining the 
latter with the requirement that all fields are charged under gauge transformations 
removes all renormalizable LV operators and the theory becomes LI in the IR within 
weak coupling. Some implications of this mechanism have been explored in [23]. 



(ii) Turn to strong coupling. Typically, the coefficients of the marginally irrelevant LV 
operators run according to 

fj. K, t t r ~ 

P 9 k lv ; , 



dlog/j, 

where g represents the coupling constant of the theory and (3 > is a numerical 
constant. If g is classically a marginal coupling it will run logarithmically (while it 



2 Indeed, any marginal symmetric tensor operator with two indices must be conserved [18, 19] and thus 
represents an alternative energy-momentum tensor. This is unlikely to happen in an interacting theory. 

3 This option is realized in the LV Standard Model Extension (SME) [20] where all dimension-3 LV 
operators turn out to be CPT odd. Therefore, imposing CPT-invariance (which is compatible with LV) 
eliminates all dimension-3 operators in the SME at once. Note that if CPT is not imposed these dim-3 
operators would seriously aggravate the fine-tuning needed to satisfy the experimental constraints on LV. 
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remains weak) and this immediately leads to the log-running for klv- To see what 
can happen at strong coupling, assume that the dynamics is such that g flows to a 
nontrivial fixed point g*. Then k lv scales like a power-law — its running gets strongly 
accelerated and the LI IR fixed point is approached at a much faster rate 4 ! In the 
context of particle physics one can then envisage the following scenario (cf. [8]): the 
theory starts non-relativistic at very high energy scale, it runs into strong coupling at 
lower energies and becomes relativistic. Finally, it enters into a confining phase at, say, 
a few TeV to give rise to the Standard Model (and, possibly, other hidden sectors) 5 . 

In this paper we consider option (ii), refining our original question as follows: What 
precise implications does the emergence of LI have for physical observables ? What kind of 
LV effects survive at low energies and how suppressed they are ? 

To have a handle on the strong dynamics, we will use the gauge/gravity correspondence. 
Namely, we are going to consider strongly coupled <i-dimensional SFT's in the limit of large N 
(number of degrees of freedom) and at large 't Hooft coupling, which are believed to be dual 
to (rf+l)-dimensional weakly coupled gravity theories defined on scale-invariant backgrounds. 
We are going to consider various deformations of these SFT's which in particular won't leave 
them strictly scale-invariant anymore, but we will continue to refer to the resulting LV field 
theories as SFT's for short. To study when and how LI can emerge in such LV field theories, 
we consider two distinct setups: 

1) The SFT is defined as a (Lorentz invariant) CFT coupled to another 'fundamental' 
Lorentz Violating sector. The two sectors are assumed to be 'weakly' coupled at some 
UV scale, i.e., with a controllably small coupling constant. This can be realized holo- 
graphically as a Randall-Sundrum type model [25] with an Anti-de-Sitter (AdS) bulk 
geometry and LV boundary conditions in the UV. Whether LI emerges (or 'persists') 
in the IR in this type of model is already nontrivial. The answer depends on the type 
of coupling between the CFT and LV sector, and we will see that LI does emerge 
provided the coupling is relevant. 

2) The SFT is by itself non-relativistic in the UV. To fix ideas, we will concentrate on 
the case when at high energies the SFT displays anisotropic scaling characterized by 
the dynamical exponent z > 1. It is conceivable that such anisotropic SFT's admit 
relevant deformations which generate a flow to a LI fixed point in the IR 6 . In these 



An alternative way to accelerate the running is to increase the number of species contributing into the 
renormalization of K^y [8]. 

5 This scenario resembles the 'walking technicolor' idea [24] invoked for the solution of the gauge hierarchy 
problem. 

6 An example exists already within free-field theory: the single-species free Lagrangian 

- 2 </>(-&)*</> 

C = ^ ~ M 2 z -2 ■ 

where A = didi is the spatial Laplacian and z > 1, exhibits anisotropic scaling. It admits a deformation by 
a relevant operator c 2 (j)A(j), which drives the theory in the IR to be LI with an emergent speed c. 
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cases, the theory exhibits an emergent LI. This can be realized holographically as a d+1 
space-time that approaches a Lifshitz geometry [26, 27, 28] in the UV. For a certain 
range of z, these fixed points do admit relevant deformations that trigger a flow to LI 
fixed point. These flows are dual to 'domain-wall' geometries that interpolate between 
Lifshitz in the UV and AdS in the IR. 

In both cases we assess the emergence of LI using two types of observables. First we 
analyze the two-point functions of scalar operators and show that they interpolate from a 
non-relativistic behavior at short distances to the LI form at large scales. Second we explore 
the dispersion relations of bound states in the confining phase of the theory, which we model 
by cutting off the geometry with an IR brane. This will allow us to study the scaling of the 
LV effects in the dispersion relations with the IR cutoff by tracing their dependence on the 
brane location. 

In this paper we do not attempt any applications of our construction to realistic model 
building. Rather, our aim is to study the basic aspects of the mechanism for emergent LI in 
strongly coupled theories. 

The paper is organized as follows. In Sec. 2 we consider the Randall-Sundrum type 
model, starting with the analysis of the correlator in Sec. 2.1 followed by the study of the 
bound state spectrum in Sec. 2.2. In Sec. 3 we turn to the model of the Lifshitz flow. 
We introduce the relevant geometries in Sec. 3.1, study the two-point function of a probe 
operator in Sec. 3.2 and the bound states in Sec. 3.3. Section 4 is devoted to discussion. 
Three appendices contain technical details. 

2 Randall — Sundrum model with Lorentz violation on 
the UV boundary 

We consider a slice of AcLS^+i-space with the metric 

ds 2 = (^j {-dt 2 + dx 2 + du 2 ) , 

where L is the AdS radius and i runs from 1 to d — 1. The two boundaries are located at 
uuv and um = A" 1 . Without loss of generality we can fix uuv = L. In this slice we put a 
probe scalar field <fi with the action 

S =k I <**"W=g( - S^Om^ - ^ 2 2 ) 

d (2) 
+ \( d 2 + c 2 (-A/ 2 )A - aly)^ f d d x nu 2 , 

A Ju=L z Ju=A- 1 

where we have used that the induced metrics on the branes are flat. Here A = didi denotes 
the spatial Laplacian, and c 2 is an arbitrary positive function. When c 2 ^ 1 the kinetic 
term on the UV brane explicitly violates the Lorentz invariance (LI). The parameter /, with 
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dimensions of length, sets the scale of this violation. It is natural to assume it to be of order 
L, though its precise value will be irrelevant for us. The parameter b will be also assumed 
to be of order L, unless stated otherwise. For simplicity, we have not included an induced 
kinetic term for on the IR brane. Still, we have included the bare mass terms on both 
branes, which will allow us to obtain massless modes by an appropriate tuning of the masses. 

Let us work out the (i-dimensional dual of the theory (2). The bulk part describes a CFT 
containing a scalar operator with dimension 7 



where 



dim 0$ = ^ + v , (3) 



j + (PL? • (4) 

The CFT is deformed both in the IR and the UV. By the standard rules, the IR deformation 
corresponds simply to the introduction of a confinement scale A. In the UV the CFT is cut 
off at the scale L -1 , where it is coupled to an elementary scalar 4>(x) with Lorentz violating 
(LV) action. The latter corresponds to the boundary value of the scalar field (f>(x, u) in the 
AdS picture. Thus the action of the dual theory has the form, 

S = Soft + l -jd d x^{- d 2 + c 2 (-A/ 2 )A - fl 2 uv )4> + j d d x 00^,, (5) 

The coefficient in front of the last term deserves an explanation. To obtain it one first 
absorbs the factor 1/6 in front of the bulk action in (2) by rescaling the field = y/b(p. 
Then the interaction of with the dual CFT does not contain any factors of b. Going 
back to the canonically normalized four-dimensional field one concludes that the coupling 
between and CFT is proportional to 1/ y/b. The dependence on L is then reconstructed on 
dimensional grounds. Note that the limit b — > corresponds to a strong interaction between 
the elementary scalar and the CFT, while b — > oo, on the contrary, gives weak coupling. 

Importantly, the interaction between the elementary scalar and the CFT can be relevant 
or irrelevant, depending on the parameter v. Indeed, the dimension of the scalar is d/2 — 1, 
implying that if < v < 1 the last term in (5) is relevant, for v — 1 it is marginal and 
for 1 < v irrelevant. Correspondingly, the physics is expected to be quite different in these 
three cases. The first case (0 < v < 1) is most interesting for our purposes. In this case the 



7 If < v < 1 the bulk theory admits an alternative interpretation as a CFT with the dimension of the 
dual operator dimO^ = d/2 — v. The value of the field on the UV brane then corresponds to the operator 
itself and the dual action (5) is replaced by 



Scft + I d d xO^{ - 8 2 t + c 2 (-Al 2 )A - f uv )0^ . 



In what follows we stick to the 'standard' interpretation based on (3) as being more general and covering 
also the case v > 1. But the analysis of Sees. 2.1, 2.2 applies without changes to the alternative picture as 
well. 
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interaction between the scalar and the CFT is relevant and thus becomes strong in the IR. 
In other words, the scalar becomes part of the strongly interacting sector. According to the 
discussion in the Introduction, one expects the theory to flow to a LI invariant point. We 
will see that this expectation is indeed confirmed. The marginal case v — 1 is more subtle. 
We will see that the theory still flows to a LI fixed point in the IR but the LV corrections 
die away only logarithmically. Finally, for 1 < v we expect that the interaction between the 
scalar and the CFT becomes less and less important when flowing down to the IR, so that 
at low energy we obtain two essentially decoupled sectors: the Lorentz invariant CFT and a 
Lorentz violating scalar. This case is studied in detail in Appendix A, where we show that 
this expectation is essentially correct, though the actual picture is somewhat more subtle. 



2.1 Propagator of the elementary scalar 

To warm up let us consider the setup with the IR brane sent to infinity which corresponds to 
the removal of the IR cutoff, A — > 0. The quantity we are interested in is the propagator of 
the elementary scalar. To obtain it we fix the value of the scalar field at the UV boundary, 
(f)(x,u = L) = (p(x), and integrate out the bulk degrees of freedom. Upon performing the 
Fourier decomposition along the (i-dimensional coordinates, 



<f>(t,x.,u) oc w ,k(w)e 



and rotating to the Euclidean signature the bulk solution reads, 

, ( v ( u\ d/2 K u (p E u) - 



Here w = —iu is the Euclidean frequency, p E = Vw 2 + k 2 and K u (z) is the Macdonald 
function. Substitution of this expression into (2) yields the quadratic boundary action, from 
which one reads off the propagator of the elementary scalar, 



G^(w,k) = (0-™-k0™,k) oc 



2 . 2,2 . 2 . PeKu-APeL^ ' 

W +Ck +fiuV+ bK„(p E L) , 



(6) 



where 



and we have used the relation, 



2 -2 \ ( d 

»uv = Huv- U {- 2 -v 



(7) 



K' v {z) = -K v - 1 (z)--K v {z). 



For high frequency and momentum, w,k ^> the last term inside the brackets in (6) is 
linear in p E and can be neglected compared to the first two terms. Thus we recover the LV 
propagator given by the brane part of the action (2). 
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On the other hand, the behavior of the propagator at low momenta is affected by the 
interaction of with the CFT. To be able to probe this regime we will assume that \i\jy 
vanishes, so that <j> describes a gap less mode 8 . The analysis is different depending on whether 
v is greater, equal or smaller that f . For v > 1 the leading term in the ratio of the Macdonald 
functions at peL <C 1 is quadratic in p E and the propagator reads 



Q~ x oc (f + x)w 2 + (c 2 + x)k 2 



where 



x = 



(8) 



(9) 



26(i/ - f) ' 

We see that, while the coefficients in front of the frequency and momentum have received a 
finite renormalization from the CFT, they are still different from each other implying that 
violation of LI persists down to low energies. This is consistent with the fact that for v > 1 
the interaction between the elementary scalar and the CFT is irrelevant. 

It is worth noting, however, that the CFT contribution grows as we increase the bare 
coupling between and the CFT (recall that the latter coupling is proportional to 6" 1 / 2 ). 
As a result the 0-propagator is driven closer to the LI form. This illustrates the general 
statement made in the Introduction that strong coupling facilitates the emergence of LI. To 
avoid confusion we stress that in this case the strong coupling is not achieved by the RG 
running, but simply by setting the bare coupling constant in the UV to be large. From the 
AdS perspective this corresponds to increasing the coefficient in front of the bulk Lagrangian, 
which in the model (2) is just a fixed parameter. 9 

We now turn to the case v < 1 when the coupling between and the CFT is relevant. 
Expanding the Macdonald functions at p E L C 1 we obtain, 



1 ocw 2 + c 2 k 2 + 2 l ' 2 ^ PEL) 



2v 



r(i 



V) 



bL 



r(iz) 



i - 



r(i 



)(p_eL\ 
r(i + i/)V 2 ; 



(10) 



The last term dominates at small momenta 10 . Remarkably, it is LI. Thus we conclude that, 
consistently with our previous expectations, the 0-propagator flows to the LI form in the 
IR. From (10) we see that the corrections to LI scale as (peL) 2 ^ v \ i.e. they are power-law 
suppressed. Note that the dimension of at the IR fixed point differs from canonical and is 
given by d/2 — v. With this renormalized value of dimension the interaction between and 
the CFT in (5) becomes exactly marginal, as it should be. 



8 In fact, it would be enough to assume [i\jv -C L -1 . Setting fiuv — is just convenient to simplify the 
formulas. 

9 One can envisage a more general setup where this coefficient is rendered dynamical and varies along 
the holographic coordinate u. This would correspond then to the RG running of the coupling which, if it 
increases towards IR, would trigger the emergence of LI. We will briefly come back to this possibility in 
Sec. 4. 

10 Note that due to the factor in the square brackets this term has a smooth limit at v — > 0. 
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The case v — 1 requires a separate treatment. In this case the propagator at small 
momenta reads, 

^W + c 2 A; 2 + ^(-log^- 7 ) , (11) 

where 7 is the Euler number. The last term still dominates and thus the IR propagator is 
LI. However, the LV corrections are only logarithmically suppressed. 



2.2 The bound state spectrum 

We now come back to the two-brane model. Introduction of the confinement scale leads 
to the appearance of the bound states, which in the holographic picture correspond to the 
Kaluza-Klein (KK) modes. Our goal is to analyze the dispersion relation of these modes. 
In this subsection we concentrate on the most interesting cases v < 1 and v — 1 where LI 
emerges at low energies automatically due to the RG flow. The case v > 1 is considered in 
Appendix A for completeness. 

The profile of a KK mode with given frequency and momentum is expressed in terms of 
the Bessel functions, 

&,,k = A lU d l 2 J v { P u) + A 2 u d / 2 Y v (pu) , 

where 11 

p = Vw 2 - k 2 (12) 
and Ai, A 2 are arbitrary constants. At the branes one obtains the boundary conditions 



= , (13a) 



u=L 



= . (13b) 



This translates into the system of linear equations for the coefficients A 1: A 2 that has non- 
trivial solutions provided its determinant vanishes, 

J uv Y IR -Y uv J IR = 0. (14) 

Here we have introduced the notations (Z stands for J or Y), 

Zvv = Z^ipL) + {UJ2 ~ ^ ~ ^ uv)h Z„{pL) , 

P 

Z« = Z„_!(p/A) + ^M^Z v (p/A) , 
where \i\jy is given by (7) and 

1 (d 



2 -2 
^IR — VlR 



bL\2 



— v 



11 Throughout this subsection we work in the Lorentzian signature. 
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The equation (14) determines the dispersion relations of the KK modes. 

In what follows we restrict to the case nuy = fj,m = that guarantees the existence of 
a gapless mode in the KK spectrum. Besides, we are interested in the modes with four- 
momenta small compared to the UV cutoff, pL <C 1, and thus we can Taylor expand the 
corresponding Bessel functions in (14). We first consider the relevant coupling case < v < 1. 
After a straightforward calculation Eq. (14) can be cast into the form, 

pY'j^MA) _ fc (LA)2( i-,) (c2 _ 1) ^ j (16) 



2A; ./^(p/A) 2L v ' v ' k 2 



where 

X-x(p/A) r(H 



r l-i/i 



r(i-z/) /pL x 2i 



If (pL) 2u <C 1, which holds automatically for v not too close 0, the second term in the square 
brackets can be dropped, so that J v -\ becomes simply proportional to J v -\. For any value 
of the spatial momentum k equation (16) implicitly determines p and hence the frequencies. 
Only the r.h.s. of this equation violates Lorentz invariance. Note that it is suppressed by 
the ratio of the IR and UV cutoffs to a positive power, (LA) 2 * 1 -"). In other words, the LV 
effects vanish as a power-law when the IR cutoff is lowered. 

To see what happens in more detail, it is useful to plot the l.h.s. of (16) as a function 
of p/A. The qualitative shape of this dependence is shown in Fig. 1. We see that it is 
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Figure 1: The l.h.s. of Eq. (16) as function of p/A (qualitative plot). 

composed of an infinite sequence of branches, the first branch starting from zero and going 
to +oo, whereas all other branches vary from — oo to +oo as p/A changes between the 
adjacent roots of J v -\. Each branch corresponds to a separate KK mode and describes its 
dispersion relation through the dependence of the effective mass of the mode p on the spatial 
momentum k. It is clear from (16) that this leads to the dispersion relations of the form, 

u 2 n = k 2 + A 2 F n ((c 2 -l)k 2 /kl v ) , (17) 
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where 

k LV = A(LAy~ l , 

and F n are dimensionless functions. Note that the scale kiv lies between the confinement 
scale A and the UV cutoff L^ 1 . 

Let us first analyze the dispersion relations of the massive modes. As long as 

|c 2 -l|A; 2 < k 2 LV , (18) 

one can expand the function F n in (17) which yields the dispersion relation as a power series 
in momentum, 

h.21 

< = A 2 (jf I) 2 + (1 + 25c n )k 2 + Ys ei Jj2T2 ■ (19) 

l>2 M n,2l 

Here ji^ v , n — 1, 2, . . ., is the nth positive root of J X - V and, assuming that LV in the UV is 
of order one, |c 2 — 1| ~ 1, the rest of the parameters are estimated as 

5c n ~ (LA)* 1 "") , M n>2l ~ A(LA)'^- 1 )/^ 1 ) ; 

€[ = ± stand for the signs of the corresponding terms. This expression deserves a few 
comments. First, we see that the deviation of the coefficient in front of the k 2 term from 
one is power-law suppressed by the small factor (LA) 2 ( 1_ ^. Second, contrary to the naive 
effective field theory logic, the mass parameters M n ^i suppressing different powers of the 
momentum are not equal. This is not surprising as these mass parameters are generated 
dynamically from two very different scales: the UV scale of Lorentz violation L~ l and the 
confinement scale A governing the physics of the bound states. While M n ^i are always higher 
than A, they can lie below or above L~ l depending on the value of v. For example, the mass 
scale M n 4 appearing in the quartic term is higher than L _1 if v < 1/2. On the other hand, 
another statement of the effective field theory still holds. The mere fact that (19) is a Taylor 
expansion implies that, as long as it is valid, the contribution proportional to k 21 is less 
important than all previous terms. In particular, within the validity of (19) deviations from 
LI are dominated by Sc n . 

When the spatial momentum reaches kiv the expansion (19) breaks down and the de- 
pendence of the frequency on k becomes complicated. However, when k further increases, F n 
stabilizes at the nearest root of J v -\ (depending on the sign of c 2 — 1 the latter value is 
either larger or smaller than 2\-v)- This implies that the dispersion relations become in this 
limit relativistic, only with the mass renormalized compared to the small-momentum regime! 
Remarkably, the deviation of the mode velocity from 1 is suppressed at all momenta 12 , 



dk 



< (LA) 2 ^ 



12 We omit the region of small momenta where the group velocity is different from 1 due to non-zero mass 
of higher KK modes. 
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Thus we conclude that for light bound states with masses small compared to the UV cutoff 
LV effects are uniformly suppressed at any value of the spatial momentum or energy. This 
is interpreted as follows: the physical scale that determines whether LV is important or not 
is neither k nor uj separately, but rather the effective particle mass p defined by (12). 

Now we consider the gapless mode. At small momenta p <C A the r.h.s. of (16) is further 
simplified by expanding the remaining Bessel functions. This yields a linear dispersion 
relation, 



k z 



1 + 



+ ... 



where dots stand for corrections of higher powers in k 2 similar to (19). This expression is 
valid under the condition (18). For k > k^y the situation depends qualitatively on whether 
c 2 is bigger or smaller than 1. In the former case the gapless mode behaves similar to the 
massive ones. Namely, its effective /c-dependent mass is real and stabilizes to p — Aj^\ at 
high momenta, see Fig. 1. Thus LI is approximately satisfied even by very energetic particles. 
On the other hand, for c 2 < 1 the effective mass found from (16) is purely imaginary. At 
k~^> kiv the dispersion relation reads, 



k 



4 

I 2 



r(i + ") ( 
r(i-i/) iv + (l - c 2 )bLk 2 



c 2 )bLk 2 



l/u 



We see that when k approaches the UV cutoff L^ 1 deviations from LI become of order one. 
It may even seem that the dispersion relation turns over and u 2 becomes negative signaling 
an instability. However, this is not true: at k > L^ 1 Eq. (16) is no longer valid because the 
condition \pL\ <d 1 is violated. One has to go back to the exact Eq. (14). A careful analysis 
shows that at k ^> L^ 1 the dispersion relation approaches the 'bare' non-relativistic form 
Uq = c 2 k 2 . We conclude that, within the considered setup, the 'superluminal' case c 2 > 1 
leads to more efficient emergence of LI than the 'subluminal' one c 2 < 1. 

For the case v — 1, corresponding to the marginal coupling between the LV scalar and 
the CFT, the analysis proceeds in essentially the same way as above. The only difference is 
the necessity to keep track of the logarithmic terms in the expansion of the Bessel functions. 
This leads to the following equation for the KK spectrum at pL <C 1, 



V 



nY (p/A) p_ _ b_ 

2J (p/A) ° g 2A 7 L 



-p 2 \ogLA = y(c 2 -l)k 2 

Lj 



The dependence of the first term on the l.h.s. on p is similar to that plotted in Fig. 1, being 
again composed of a series of branches of tg-like shape. The second term involving the large 
logarithm — log LA shifts every branch upward. This implies the following picture for the 
massive modes. At small k their dispersion relations start from 



m. 



A^+OdlogLAr 1 )) 



If the LV scalar is 'subluminal', c 2 < 1, the /c-dependent mass p(k) always stays close to this 
value, approaching Aj^ at k — > oo. On the other hand, in the 'superluminal' case c 2 > 1, 
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the effective mass goes to the next root of the Bessel function, p(k) — > Aj'q at k — >■ oo. 
The transition occurrs at (c 2 — l)k 2 ~ k\ v with 

k LV = A^-log LA . 

This is still parametrically higher than A, but the hierarchy is only logarithmic. The devia- 
tions of the mode velocities from 1 are suppressed by the same logarithm. 

For the gapless mode one again starts from the regime p < A where one obtains linear 
dispersion relation, 

(20) 

At A — > the velocity of the mode tends to 1 implying emergence of LI. However, the 
approach to the LI regime is only logarithmic and thus too slow to be relevant for realistic 
particle physics models. The asymptotics of the dispersion relation at high momenta !;> kiv 
depend on the sign of c 2 — 1, the analysis being completely analogous to the case of relevant 
coupling. If c 2 > 1 the mode becomes effectively massive with constant mass while for 

c 2 < 1 the dispersion relation approaches uj 2 = c 2 k 2 . 

Let us attempt an intuitive understanding of the surprising result of this section that 
the bound states display the emergent relativistic dispersion relation even at large momenta 
k ^> kiv- Actually, we have found this behavior for all KK modes except the zero mode 
in the subluminal case c < 1. In the 5D picture, the reason behind this observation is that 
the purely-LV term on the UV brane (c 2 — l)k 2 (p 2 5(u — uuv) ac ts as a k— dependent mass 
localized on the brane. As a result, at large k and c > 1 the KK wavefunctions are repelled 
from the brane and the LV effects vanish. For c < 1 the effective squared mass is negative 
and the zero-mode is attracted towards the UV brane leading to the strengthening of LV. 
On the other hand, the excited KK modes, being orthogonal to the zero-mode, are still 
repelled from the brane and their k— dependent effective mass saturates in the UV at the 
value corresponding to the Dirichlet boundary conditions. 

This has the following CFT interpretation. The emergence of LI from the RG flow means 
that the theory possesses a preferred frame which can only be 'noticed' by short- distance 
probes. The excited KK modes map to resonances — the states which are composite at all 
momenta and have a finite size of the order of the confinement scale A. Thus, they simply 
cannot probe the preferred frame. The zero-mode, instead, can be viewed as a mixture of an 
elementary and a composite state. The size of such an object can be momentum-dependent. 
For c < 1 the size shrinks with k and the state becomes mostly elementary. It is now able 
to probe the UV physics and its velocity approaches the UV value c. For c > 1, on the 
contrary, the admixture of the elementary component never becomes large and the state 
remains mostly composite. Then just like for the resonances the emergent relativity persists 
even for large k. 



2 i 2 

Wq = k 



1 + 



(Lib) log (LA) 
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3 Lifshitz flows 



While the model considered in the previous section provides a particularly simple realization 
of the dynamical emergence of LI through strong coupling, one may object that the LV sector 
on the UV brane is rather ad hoc. Admittedly, it is not straightforward to couple the action 
of the UV brane in (2) to the metric, as required for the fully dynamical formulation 13 . In 
this section we study another model that does not suffer from this drawback, being generally 
covariant by construction. It was introduced in [26, 27] and allows to describe RG flows 
between non-relativistic fixed points with anisotropic (Lifshitz) scaling in the UV and LI 
field theories in the infrared. After reviewing the model and the relevant geometries in 
Sec. 3.1, we address the physical aspects of emergent LI in Sees. 3.2, 3.3. 



3.1 The flow geometry 

This subsection closely follows the presentation of [28]. Consider a (d+l)-dimensional theory 
containing gravity with negative cosmological constant and a free massive vector field 14 . The 
action reads, 



1 



16nG 



d d+1 xJ= 



g 



2 



R — 2A C — -^FlnF ln —AnA n 



(21) 



where Tln = QlAn — OnAl, as usual. We look for solutions of this theory that are static 
and invariant under translations and rotations in (d — 1) spatial directions. This leads to 
the Ansatz, 15 



2 d-1 

ds 2 



(£) ( - + t, <*°? + 9 V) • A, = jLfj , (22) 

where /, g, j are functions of the coordinate u to be determined. It is convenient to introduce 
one more function h(u) parameterizing the field strength, 

J^ut — —h . 



u 2 



13 It is possible that such formulation can be achieved by introducing additional degrees of freedom that 
would define a preferred frame on the UV brane, similar to the approach of [35, 36]. 

14 The original model of [26], formulated in d = 3, includes masslcss 1- and 2-form fields with a Chern- 
Simons-type mixing. These are equivalent to a massive vector on-shell. 

15 Thc imposed symmetries also allow a non-zero value of the A u component. However, this vanishes in 
the metric (22) due to the equations of motion for the vector field. 
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Then the Einstein-Maxwell equations reduce to the set of first-order ordinary differential 
equations, 

uti = -MLgj + (d- l)h (23a) 
uf = -MLgh - d -^j -j^ih 2 - f + KL 2 ) , (23b) 

^l = -^-(h 2 + j 2 + A c L 2 )-^ (23c) 

Y = -^(h*-? + KJ?) + \. (23d) 

Note that the first three equations form a closed system for the functions g, h, j. 

Let us investigate the fixed points of this system. These are characterized by constant 
values of g, h, j, and a power-law behavior of the function /, 

/ = fou 1 ' 2 . (24) 

To understand the physical meaning of these solutions, note that the corresponding geome- 
tries possess an isometry, 

1 1-)- \ z t , x I—)- Ax , (25a) 
uh> Am, (25b) 

where A is an arbitrary number. If z — 1 one recovers the AdS space-time that by the 
standard AdS/CFT dictionary is dual to a relativistic conformal field theory. By the ex- 
trapolation of this logic, the geometries invariant under (25) with z different from 1 have 
been conjectured [26] to be dual to strongly coupled rf-dimensional non-relativistic theories 
invariant under the anisotropic scaling transformations (25a). This kind of scale invariance 
was considered for the first time in the seminal works by Lifshitz [37], hence we will refer to 
the geometries with the isometries (25) as 'Lifshitz space-times'. 

Setting the l.h.s. of Eqs. (23) to zero one observes that AdS is always a solution, provided 
the vector field vanishes, h = j = 0. By an appropriate rescaling of coordinates / and g can 
be set to 1, which fixes the relation between the AdS radius and the cosmo logical constant, 

2 

Clearly, this fixed point is Lorent invariant. Besides, for the vector mass in the range, 

d-i<( ML y<^t^l, (26) 

the system (23) has two fixed points with Lifshitz scaling, 
d-l 



2 

9± = 



2(ML)' 



2 \/A(A V\2 \ 2 



(27a) 



2 _ d(d-l) (ML) 2 (d-l) 2 . 2 _ (ML) 2 d-l (ML) 2 g 2 ± 
l± ~ 2 2 2g\ ' 3± ~ 2 2g\ ' Z± ~ d-l 
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For (ML) 2 < d - 1 only the "+" branch survives (j 2 _ determined from (27b) becomes 
negative) and for (ML) 2 > d(d - I) 2 /(3d - 4) there are no Lifshitz solutions at all. Note 
that on the "— " branch the critical exponent is confined to the interval 



1 < z < d- 1 , 



(28) 



while on the "+" branch it always exceeds 16 d—1. As we are going to explain shortly, only 
the "— " branch is relevant for our purposes. Thus we will restrict to the choice (26) in what 
follows. 

Our aim is to study RG flows connecting non-relativistic SFTs in the UV to relativistic 
CFTs in in the IR. These correspond to the solutions of the system (23) that interpolate 
between Lifshitz space-time at small values of the coordinate u and AdS at u — )■ oo. To work 
out if such 'domain wall' solutions can be realized in the model at hand, let us analyze the 
stability properties of the fixed points. To this end we add small perturbations Sf(u), etc., on 
top of the fixed point solutions and expand Eqs. (23) to linear order in these perturbations. 

We start with the AdS fixed point. At linear order the perturbations of the metric and 
the vector field decouple. Consider first the metric perturbations keeping the vector at zero, 
h = j = 0. One finds two modes, 5f = const , 8g — and Sf, 5g oc u d . The physical 
interpretation of these modes is known: the first corresponds to the source and the second to 
the expectation value of the energy-momentum tensor of the dual CFT. Indeed, the growing 
mode represents the first term in the expansion of the AdS-Schwarzschild metric. Thus, the 
geometry that is obtained by the excitation of this mode is that of an AdS black hole. In the 
dual picture this corresponds to setting the CFT at finite temperature. As in the present 
paper we are interested in the vacuum solutions, we shall set this mode to zero. On the 
other hand, the constant mode can be eliminated by the rescaling of the time-coordinate. 

Another set of modes involve perturbations of the vector field. These behave as the 
power-law, 



To have a domain wall solution that asymptotically approaches AdS at large u we need at 
least one decaying mode. In other words, the vector field must be heavy enough, so that 



Note that this condition coincides with the lower bound in (26). It matches nicely with the 
holographic interpretation of the considered system. Indeed, the presence of the vector field 
Am in the bulk theory implies that the dual LI CFT contains a vector operator O^, whose 
dimension is given by the exponent of the growing mode, 



16 For a given value of ML in the range (26) the critical exponents on the two branches are related by a 
simple formula = (d — l) 2 . 



h , j oc u 



-a± 




a = a + > . 



(29) 



dim = d + a . 



(30) 
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The domain wall solution we are looking for corresponds to perturbing the CFT with the 
time-component of this operator (and in this way breaking LI). For the LI fixed point to be 
IR attractive, the operator must be irrelevant, dimO^ > d, which is equivalent to (29). 
Note, however, that the upper limit on the allowed masses of the vector field in (26) implies 
that a is bounded from above and thus the dimension of cannot deviate parametrically 
from d. For example, for d = 3,4 we have that a is smaller than 0.13 and 0.35 respectively. 
In fact, a is less than 1 for all d < 6. 

Let us turn to the perturbations around the Lifshitz fixed point. The analysis us identical 
for the "+" and "— " branches. Linearizing the first three equations of (23) one finds three 
modes, 

Sh, Sj, 5g oc u 13 , 

where j3 can take the values, 

/3 1 = z + d-l, (31) 
£± = \ [z + d - 1 ± y/(z + d- l) 2 + 8(z-l)(z-d+ 1)] . (32) 

Inserting these into the linearized Eq. (23d) we obtain the corresponding perturbations Sf 
that we define as, 

f = ftiu 1 -'(l + Sf). 

At this stage we find a fourth mode Sg = Sh = Sj = 0, Sf = const. In [28] this mode 
together with the mode corresponding to fli (31) were interpreted respectively as the source 
and expectation value of the energy density of the system. Indeed, the energy density 
operator must be marginal, i.e. its scaling dimension must be (z + d — 1) — the inverse of 
the scaling dimension of the integration measure dtd d ~ 1 x in the action 17 — which matches 
with f3\. 

Consider now the pair of modes (32). Following the usual holographic vocabulary one 
interprets the mode corresponding to /3_ as the source of a certain operator Oj^ in the 
Lifshitz theory, and the /3 + -mode as its expectation value [28]. Consequently, (3 + sets the 
scaling dimension of the operator. For 1 < z < d — lwe have j3 + < z + d — 1 implying that 
the operator is relevant, for z = d — 1 it is marginal, while for z > d — 1 it is irrelevant. 
Thus the Lifshitz fixed point can be UV attractive only in the range (28). This excludes 
from consideration the "+" branch of fixed points. The geometrical counterpart of this 
statement is that if and only if 1 < z < d — 1 both exponents (32) are positive implying 
that all perturbations decay towards the asymptotic UV boundary u — >■ and the Lifshitz 
space-time is UV attractive 18 . 

17 In our conventions the scaling dimension of the spatial coordinates is —1. 

18 The situation is somewhat trickier for the marginal case z = d — 1. The mode corresponding to /3_ 
is then constant at the linear level and one has to take into account the non-linear corrections to lift this 
degeneracy [30]. One finds that due to the non-linearities the mode is driven to zero logarithmically at 

u — y 0, 

Sh, Sj, Sg oc | logu| _1 . 
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Figure 2: The vector field and metric profiles for the domain wall solutions interpolating 
between the Lifshitz and AdS spacetimes. The values of the parameters are d — 4 and 
z = 1.5, 2, 2.5, 3 (from bottom to top). One may notice from the graphs that the approach 
to the UV asymptotics is slower for z = 3 than for the other cases, signaling the presence of 
the marginally relevant mode (see footnote 18). For other values of d the curves are similar. 



This suggests a simple strategy to search for domain wall solutions with desired properties. 
Start at large u from AdS slightly perturbed by the mode oc u~ a and integrate Eqs. (23) 
towards the decreasing value of the coordinate. As the Lifshitz fixed point on the "— " branch 
is absolutely stable at u — > 0, one may expect that the solution will be attracted to it at 
small u. This is indeed confirmed by direct numerical integration of the system (23), the 
resulting solutions are plotted in Fig. 2. It is worth stressing that from the dual viewpoint 
the existence of these solutions is highly non-trivial. There one starts from a fundamental 
Lifshitz SFT and adds to it the (marginally) relevant deformation Oj^. It turns out that the 
resulting RG flow is driven to LI in the IR 19 . 

Let us fix a few notations that will be used in the subsequent sections. First, from Fig. 2 
we see that the transition from Lifshitz space-time to AdS is quite sharp, which allows to 
introduce the notion of the domain wall position (this can be defined, e.g., as the value of 
the li-coordinate where the function j(u) is half of its value at u — > 0). We will denote this 
position by A" 1 . Clearly, there are domain walls with any possible value of A~ x that are 
related to each other by a rescaling of the coordinates. From the physical viewpoint, A* sets 
the scale of the transition from non-relativistic to LI regime. To simplify the formulas we 

However, this introduces an additional logarithmic scaling into the function / on top of the power-law 
dependence, 

/ oc u 1 " 2 | logu| _2 ~ , 

so that the UV asymptotic is not, strictly speaking, Lifshitz. Still, it has been argued in [31] that this can 
be interpreted as just a mild violation of the Lifshitz scaling in the dual theory by the marginally relevant 
deformation, so that for large spans in energy the physics is well approximated by the Lifshitz fixed point. 

19 Interestingly, for the opposite sign of the relevant deformation, one obtains a domain wall solution 
connecting to the "+" -branch of Lifshitz fixed points in the IR [28]. 
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will work in the units such that 

A* = 1 . (33) 

Second, we will need the coefficients of several leading terms in the expansion of the domain 
wall configuration at u — > oo. We write, 

h « hooU' 01 , j « jooU~ a , (34a) 
f^l + Uu- 2 *, g^l + 9oD u- 2a . (34b) 

In the expressions in the second line we have used that the perturbations of the metric are 
induced only at the second order in h,j as is clear from (23c), (23d). The coefficients in (34) 
are not independent. They can be determined from (23) up to a single constant. However, 
the corresponding relations are rather cumbersome and we are not going to use them in 
our analysis. The overall normalization of the coefficients in (34) is connected, in turn, to 
the domain wall position. With the convention (33) all these coefficients are naturally of 
order one. Finally, the numerical solution shows that all functions f(u), g(u), h(u), j(u) are 
monotonically decreasing at all u > 0. In particular, this implies that the coefficients /oo, 
fi'oo, ^oo, Joe are positive. 

A comment is in order. The Lifshitz space-times and the domain wall solutions described 
in this section are known to be geodesically incomplete at u — > oo where they suffer from 
a specific type of singularity [32, 33]. Despite the fact that all scalar curvature invariants 
are bounded, the tidal forces in a geodesically moving frame diverge at u = oo in a finite 
proper time. We will not attempt to provide any resolution or physical interpretation of this 
singularity 20 . Instead, we think about the simple setup based on the Einstein-Proca action 
(21) as a convenient toy model illustrating general features of RG flows with emergent LI, 
which are expected to be valid in more sophisticated backgrounds with holographic Lifshitz 
duals (e.g. like the supergravity solutions of [28]). 



3.2 Probing the flow with the two-point correlator 



To study the physical properties of the holographic RG flow we make use of a probe scalar 
field that we add to the bulk action. This implies that the dual strongly coupled theory 
contains a scalar operator 0$. The coupling of <fi to gravity is taken to be minimal for 
simplicity. On the other hand, we introduce a non-minimal interaction with the Lorentz 
breaking vector Am in order to capture the species-dependent effects of LV. This leads to 
the action, 



= - I d d+1 Xs 



(MLf 



■A M A N 



d M (j>d N (fr - fi 2 (p 2 



(35) 



Due to the ^-coupling the effective metric felt by the field is different from g MN , so that the 
0-modes propagate with a speed different from 1. For £ > they are "subluminal" (speed 
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This issue was discussed in [34]. 
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smaller than 1), while for 21 £ < they are "superluminal" (speed greater than 1). Note that 
we could add several probe fields with different values of the parameter £ meaning that this 
parameter encodes the species-dependent information. 

Using the operator one can construct various observables probing the properties of the 
RG flow dual to the domain wall solution. In this section we concentrate on the two-point 
function, 

G^u, k) = (O^cu, k)O (-w, -k)> . (36) 

At small distances, where the relevant perturbation of the Lifshitz fixed point is not effective, 
one expects this correlator to obey the Lifshitz scaling. On the other hand, at large distances, 
i.e. at small energies and momenta, its behavior must be dominated by the LI infrared fixed 
point and therefore it must exhibit approximate LI. Our aim will be to verify the latter 
assertion. Note that it is not immediately obvious, given that the holographic prescription 
for the calculation of the correlator involves taking the limit of certain quantities at u — >■ 0, 
i.e. deep inside the Lifshitz region. Additionally, we will find the form of the LV corrections 
to the correlator at low energies. 

The calculation is based on the solutions of the equation of motion for the scalar field 
obtained from (35). This reads in the Fourier representation, 



i + ( 1 ---—^-g {-y-w +k + — H0 = O, (37) 



where primes denotes the derivatives with respect to u and we performed the Wick rotation 
to the Euclidean frequency, w = —iui. This equation cannot be solved exactly because we 
do not have analytic expressions for the functions f,g,j of the domain wall background. 
However, for our purposes the explicit form of the 0-modes will not be required. It suffices 
to know their qualitative features, which we now discuss. 

In the Lifshitz region m < 1 the background functions can be replaced by their asymptotic 
form 22 (24), (27), where the "-" branch must be taken 23 . Then one can choose two linearly 
independent solutions of (37) that have the form of a simple power-law at u — >■ 0, 



V± = u 



x± , u — y , (38a) 



d + z- 1 / fd + z-V 1 



\± = — n ±\ [—^ +{9^Lf. (38b) 



21 Thc absolute value of £ must be small enough for the combination 1 + £,j 2 (u) to be positive on the 
domain wall solution. Physically, this requirement means that the term with time-derivatives in the action 
(35) is everywhere positive, so that the spectrum of excitations does not contain ghosts. 

22 We assume that z is strictly smaller than d — 1 to avoid the complications appearing in the case of 
marginally relevant deformation. 

23 The expressions in terms of the critical exponent are, 



lz 2 + (d-2)z + (d-l) 2 . / d{d-lf{z-l) 



d(d-l) ' J ~ V z 2 + (d-2)z+{d-l) 2 ' 
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Here the subscript "a" reflects the fact that these solutions are defined by their behavior in 
the Lifshitz part of the space-time, to which we will refer as the "a-region" for short, see 
Fig. 3. 

We need a linear combination of the solutions (38) that decays at infinity, u — > oo. This 
will be denoted by (pb(u), being defined by the behavior in the "6-region" u ^> 1. We write, 

Uu) = T_cf> a ,_(u) + T + a , + ( M ) , (39) 

where the coefficients T T depend on the frequency w and momentum k. According to the 
dictionary of the holographic correspondence the amplitude T_ is identified with the source 
for the operator O^. The correlator (36) is then expressed through the variation of the 
classical action (35) evaluated on the solution (39), 

5 2 S 

^ (w ' k)oc ^Kk)^(-^-k) +c - t - 
iff 

If u-*o \ g 

The proportionality signs here emphasize that we do not keep track of the overall normal- 
ization of the correlator focusing only on its frequency and momentum dependence. The 
notation "ct." stands for the counterterms that must be included to remove the divergen- 
cies that arise in taking the limit indicated in the second line. Using (39) we obtain, 



G^W, k) OC lim [U 1 - d -<t> a ,-(t>'a - + C.t. 
u— s-0 \ g ' 



+ (A + + A-)/o 

9- T-(w,k) 



Only the first term diverges at u — > and must be cancelled by the counterterms. Impor- 
tantly, this term is polynomial in w and k, meaning that it is local in space and time. Indeed, 
at u — y the contributions containing w, k in Eq. (37) are subleading and the solution <p a _ 
can be obtained by iterating them in a perturbative expansion. As a result <fi a ~ is obtained 
in the form of a power series in w 2 and k 2 . Only a finite number of terms in the series lead 
to divergencies when substituted into (40), implying the above assertion. Thus we conclude 
that the first term in (40) does not affect the correlator of operators taken at different space- 
time points, which is given by the non-analytic part of G</>(w,k). In what follows we will 
discard polynomial contributions into the correlator. Thus, we are left with the second term 
in (40). In other words, the correlator is given by the ratio of the coefficients in the linear 
relation (39). The situation is exactly the same as in the standard case of the relativistic 
AdS/CFT correspondence, cf. [38]. 

Before concentrating on the case of small frequencies and momenta, which is of the main 
interest to us, we briefly discuss the opposite regime w,k ^> 1. The decaying solution fa is 
predominantly localized at u < 1/k and thus for large momenta it is entirely controlled by 
the Lifshitz region m < 1. By an appropriate rescaling of variables one finds that it has the 
form, 

i + a 2 , w 2 

fh 2 -' 2 k 2 > 



% I 9-ku ; 2z _ 2 — , g~fiL 
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Figure 3: The regions on the w-axis that enter into the construction of the solutions to 
Eq. (37). The case of small frequency and momentum is shown. 

where rji, 772) is the decaying solution of the equation 

4"+ 2 ~ d ~ z $ - Ux 2z - 2 + 1 + 4V = ° • 

X \ X 2 J 

Substitution into (39), (40) yields, 

For the case d — 3, z — 2 the function Q was found explicitly in [26]. We see that at large 
momenta the correlator obeys the Lifshitz scaling with the critical exponent z. The overall 
scaling of the correlator with momentum implies that the dimension of the operator in 
the Lifshitz fixed point is A+ (recall that A + — A_ = 2A+ — z — d + 1), which is consistent 
with the holographic interpretation of the solution (see (38)) as the expectation value 
of Otj,. Note also that the dependence of the correlator on frequency and momentum is not 
universal: it contains the species-dependent coupling £. 

We now study in detail the behavior of in the low-energy regime w, k <C 1. In this 
case the structure of the solution 4> b is more complicated: it feels all of the domain wall 
configuration and has a tail extending deep inside the 6-region. In the latter region the 
space-time is approximately AdS, so one can find 0^ analytically by expanding in the small 
deviations of the background functions /, g, h from their AdS form. On the other hand, at 
u <C w~ l ,k~ l (we will call this "a'-region") the terms proportional to the frequency and 
momentum in Eq. (37) can be treated as perturbations and the solution can be expanded as 
a power-series in w 2 , k 2 . This suggests the following perturbative strategy to construct 0&: 
find the two expansions above to a given order in the corresponding domains and match them 
in the overlap 1 u it; -1 ,/? -1 . The latter overlap will be referred to as the "c-region". 
The regions involved in the analysis are depicted in Fig. 3. 

Consider first the leading order. In the b- region one sets f = g = l, j = and obtains, 

0[ o) = u d l 2 K v {p E u) , (41) 
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where we use the same notations as in Sec. 2.1 and v is defined in (4). In the c-region this 
function is matched to the linear combination, 

^\u) = U^^l(u) + U^^l(u), (42) 

where <j)^± are the solutions of Eq. (37) with w and k set to zero. They are fixed by their 
asymptotics at large u, 

<$l « u d l 2±u at u > 1 . (43) 
Using the expansion of the Macdonald function at small values of the argument yields, 

Jj(0) = 2 Tu-l T{Tv)p ±, (44) 

To find the coefficients T± entering into (39) we need to rotate from the basis 0^,°±(' u ) i n 
the space of solutions to Eq. (37) with w = k = to the basis <fi^± defined by the asymp- 
totics (38). We write 

^ = ym+m + ym^m , ^<»> = v m^_ + y^m . (45) 

The coefficients V±\_ are sensitive to the shape of the domain wall background and cannot 
be found explicitly. For our purposes it is sufficient to know that they do not contain any 
dependence on w or k. Combining the linear relations (45), (42) and evaluating the ratio of 
the resulting coefficients 7± we obtain, 

r(-v) |\/ (0) l 

Q+(w, k) oc 2- 2 '^^_J_ . $ + . . . , (46 ) 

where | | is the determinant of the 2 x 2-matrix composed of the coefficients 24 V^°]_ and 
dots stand for the terms with higher powers of p 2 ^ that are suppressed at pe <S 1. In deriving 
(46) we have subtracted a constant (frequency and momentum independent) contribution. 
The expression (46) is manifestly Lorentz invariant. It has the right dependence on the 
momentum for the two-point function of the scalar operator with dimension (3), which is 
indeed the dimension of 0$ at the relativistic fixed point. Thus we have checked that at 
large distances the behavior of the correlator changes from Lifshitz to LI, as expected. 

Next we want to investigate the LV corrections to (46). This requires going to the next 
order in the perturbative scheme for the solution 0{>. The calculation is rather tedious and 
we present it in Appendix B. The final result has the form 25 , 

Q4W, k) OC 2- 2 - i > i_J- . ft . [1 + w 2 p -2 + 2a + ^ + ^ + ] > (4?) 



) {v^ly 



24 0ne can show that |V^| > 0, see Appendix B. This ensures that the correlator has correct analytic 
properties: its imaginary part is positive on the upper edge of the cut at p 2 = —p 2 E > 0. 

25 Recall that a has been defined in Sec. 3.1 and is related to the asymptotics of the domain wall background, 
see (34). 
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where the coefficients 61,2,3 are given in (88) and dots stand for terms that are either of 
higher power in frequency and momentum, or are LI, i.e. depend on w, k only through the 
combination p E . Note that we have kept the C 3 -term in spite of the fact that it is subleading 
compared to the (^-contribution; the reason for this will become clear shortly. One observes 
that the LV corrections are in general different for different species: the formulas (88) for 
the coefficients depend both on the parameter v related to the dimension of the operator 
and on the species-specific coupling £. 

As long as a < 1 the Ci-term in (47) dominates meaning that the scaling of the LV 
correction with energy is determined by the dimension of the vector operator O^, see Eq. (30). 
If a > 1 the CVterm becomes dominant and the LV correction scales as w 2 . Finally, for 
a > 1 + v the Ci-term is small compared to the (^-contribution and the leading corrections 
are completely controlled by the properties of the operator itself. Note, however, that 
due to the upper bound on a discussed in Sec. 3.1 this situation cannot be realized within 
the considered setup for spacetimes of low dimension. 

The structure of the corrections (47) admits nice holographic interpretation. They corre- 
spond to an irrelevant perturbation of the relativistic CFT emerging at the IR fixed point by 
the leading LV operators. Apart from the operator 0\ discussed in Sec. 3.1, we have at our 
disposal the time derivative of O^. Omitting dimensionless coefficients we write schematically 
for the perturbation Lagrangian, 

= f + aS + (°* + ^)' j+ "" <«) 

where we have restored the LV scale A* in appropriate powers. Dots stand for operators 
of higher dimension that are further suppressed by A*. Note that we have included the 
interaction of with the source J that also receives LV contributions due to the RG flow. 
Now, various terms in (47) are identified with the contribution of (48) into the perturbative 
calculation of the O^-propagator around the IR fixed point. The Ci-term comes from the 
double insertion 26 of the first term of (48) into the propagator, the C^-term is identified with 
the insertion of the second term, and the C2-term comes from the modification of the coupling 
to the source. An immediate consequence of this identification is that the LV correction of 
the CVtype is not universal in the sense that it will not appear in the observables which 
do not involve external sources. Indeed, we will see in the next subsection that this type 
of corrections does not show up in the dispersion relations of the bound state present in a 
confining version of the theory. On the contrary, the LV corrections of the Ci, C3- types are 
expected to be universal as they are intrinsic to the theory itself. 

26 A 

single insertion gives zero due to vanishing of the three-point function, (O^O^O 1 ^) = 0, in the un- 
perturbed CFT. The latter property follows from the invariance of the bulk theory under the flip of sign 
Am >-> Am- 
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3.3 Bound states 



To make contact with particle physics, the model of the previous subsection must be supple- 
mented by a mechanism providing a discrete particle spectrum. This is achieved by cutting 
the space-time with a brane located at u = A" 1 and keeping only the portion < u < A -1 . 
Note that in this in way one shields the curvature singularity at u — > oo mentioned at the 
end of Sec. 3.1. In Appendix C we show that the brane energy-momentum tensor required 
for consistent matching with the domain wall solution is a sum of a constant negative tension 
plus a contribution that satisfies the null energy condition and thus can be provided by some 
regular matter. This setup gives rise to a discrete spectrum of KK modes of the scalar field 
(f) that by the holographic correspondence are associated to bound states in the dual theory. 
From the dual viewpoint the bound states appear due to the introduction of an IR cutoff 
(confining scale) A set by the brane, to which we will refer as "IR brane" in what follows. 
We will assume the confining scale to be much lower than the scale of LV, A < A„, which 
under the convention (33) is equivalent to 

A < 1 . 

Our aim is to determine the dispersion relations of the bound states (equivalently, KK 
modes). Similarly to the model of Sec. 2, we introduce a mass /j,tr for the field <fi on the IR 
brane that will be tuned to make the lightest KK mode gapless. 

It is convenient to put the eigenmode equation (Eq. (37) with w 2 replaced by — u 2 ) into 
the form of a quantum mechanical eigenvalue problem, 

2 k 2 )^ 



(H {0) + V(u))cf> 
where the operator is defined as, 

p r d 2 (r 



H, 



(o) 



du 2 



and 



V(u) 



f 
f 2 



9 



(49) 



d-l\ d 



du 



u 



+ g- 



21 



1 



k 



Equation (49) is supplemented by the boundary conditions, 

-> at u — y , 



0' 



. 



(50) 

(51a) 
(51b) 



u=A- 



Note that we impose the Dirichlet condition in the UV. The IR boundary condition (51b) 
is due to the mass term on the brane, cf. (13b). It is straightforward to check that the 
operator is Hermitian in the space of functions satisfying these boundary conditions 
with the scalar product 



(01 1 02 



duu 1 -* g(1 + ^' 2) <j>l(u)<h(u) . 



(52) 
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The key idea of the calculation below is to treat the potential (50) as a small perturbation. 
This is certainly true for low-lying KK modes with small enough momenta. Indeed, V(u) is 
localized at u < 1, while the low-lying modes are expected to spread over the whole interval 
u < A -1 . Another crucial observation is that the leading operator does not contain any 
dependence on the spatial momentum k. Thus neglecting V we obtain in the leading order 
a set of LI dispersion relation 

2 i 2 i 2 

u n = k +m n) 

where the masses m 2 , n = 0,l,2,..., are given by the eigenvalues of -ff(o)- The LV corrections 
are then found using the standard techniques of quantum-mechanical perturbation theory. 



3.3.1 Unperturbed eigenfunctions 

To proceed along these lines we will need the matrix elements of V(u) between the eigen- 
functions of the unperturbed Hamiltonian H(o)- Thus, our first task is to find these eigen- 
functions. We focus on the light KK modes with the masses of order the confinement scale 
A. This implies m n <C 1, and we can use the same technique as in Sec. 3.2: find the solutions 
separately in the a'- region u <C l/m n , in the b- region u ^> 1 and match them in the overlap 
1 < u < l/m n (c- region). In the a'- region the mass can be neglected altogether and one is 
left with the equation H(p)(f) = 0. This equation already appeared in Sec. 3.2 and its solution 
vanishing at u — > has been denoted by Thus we write, 

4>n{u) = a n (f)^ + (u) , u < l/m n , (53) 

where the normalization coefficient a n will be fixed below. On the other hand, in the 6-region 
the functions /, g, j can be set to their asymptotic values, / = g — 1, j — 0, and the solution 
reads, 

(f) n (u) = b n u d/2 J u {m n u) + b n u d/2 Y u (m n u) , u > 1 . (54) 
Using the expansion of the Bessel functions in the c-region and the asymptotics, 

T/(0) v (o) 
AO) ^ ..d/2+u V +- ..d/2-u ^ I /crc^ 

^ + ^]yW\ u ~]yW\ u ' M>>1 ' (55) 

obtained by inverting Eqs. (45), we find the relations, 

y(0) _ y(P) 

bn = a «yj7^ 2 " r ( 1 + z/ ) m r i " ' h n = a «y^(5)| 2 ^smTTiTXl - v)m v n . (56) 

One observes that the coefficient b n is suppressed by m 2 ^ 1 compared to b n . Thus, the 
second term in (54) can be neglected when m n u > 1. As the KK masses are expected to 
be of order 27 A, this implies that the second term in (54) is irrelevant over most of the 



27 



Except, probably, the zero mode that must be treated separately, see below. 
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interval u < A 1 . Taking this into account we obtain from the boundary condition (51b) the 
eigenvalue equation, 

J u+1 (m n /k) - — [bLjj 2 IR + ^ + v J J u (m n /K) = . (57) 

At this stage, to somewhat simplify the subsequent formulas, we will impose the require- 
ment that H(q) has a massless mode. This has the form (53) in the whole interval. Using 
the asymptotic expression (55), where we neglect the second term, we find that it satisfies 
the boundary condition (51b) if 



\Ar = - I o + v 



d \ 1 
2 +V )bL 



Note that according to Eq. (57) the masses of other modes are then given simply by the 
zeros of J u+ ±. 

It remains to fix the normalization of the modes. Using the scalar product (52) we impose 



Jo J 



The integral is saturated at u ^> 1, so we use Eqs. (54) and (55) for the massive and massless 
modes respectively (where in both expressions we neglect the second term). This yields, 

a = ^2(1 + v) ^ 1+uL ^) ' b n = V2A\J v (m n /A)\- 1 , n>\. (58) 

In deriving the second expression we have used that J u+ i{m n / 'A) = 0. Equations (53) - 
(56), (58) completely specify the eigenmodes of H( ). 



3.3.2 First-order LV corrections 

We are now ready to evaluate the LV corrections to the dispersion relations of the modes. 
The first-order correction has the form, 

5< (1) = (0 n |U|0 n ) = A; 2 ^ A cW-^(V-i±^(0 n ) 2 . (59) 

The two factors in the integrand are peaked at different values of u. The combination in 
the brackets deviates from zero at u < 1, with a tail that extends to larger u falling down 
as u~ 2a (see (34)). On the other hand, the square of the wavefunction has maximum at 
u > l/m n ^> 1 and decreases as u d+2u towards smaller u. Depending on which factor 
dominates, the integral is saturated either in the IR (i.e. u ^> 1) or in the UV {u < 1). Let 
us consider the corresponding options case by case. 
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(i) If a < 1 + u the integral is saturated in the IR 28 . In this region we use the asymptotics 
(34) and the expression (54) for the wavefunction. Expanding the integrand to the leading 
order in u~ 2a we obtain, 



5^1 (1) = k\2U ~ Hl) 1 \ + lf V _ 0l A 2 ° , (60a 

2A 2 m~ 2+2a r mn/A 
(J„(m n /A)) 

The combination multiplying k 2 on the r.h.s. has the physical meaning of (twice) the cor- 
rection to the propagation velocity of the mode. Notice that it depends on the non-minimal 
coupling £ of the scalar field to the vector, as well as on the mass of the mode, implying 
that it is non-universal. At the same time it is clearly suppressed by the RG flow. Indeed, 
assuming m n ~ A and restoring explicitly the LV scale A*, we see that all corrections to the 
velocity are of order (A/A*) 2 ". 

One may wonder if the suppression persists for heavier modes. To answer this question 
let us consider the limit of Eq. (60b) for m n ^> A (but still assuming m n < 1). Depending 
on the value of a the integral entering this expression converges or not when the upper limit 
tends to infinity; this gives two cases, 

^«*'W.-£&)x{^ i ll 1 '] (6!) 

where we have omitted irrelevant numerical factors of order one. We see that in the first case 
the correction to the mode velocity does not scale with the mass, while in the second case 
it grows. However, in both cases the corrections are small as long as m n does not exceed 1, 
where the formula (61) stops being applicable anyway. 

(ii) For a > 1 + v the integral in (59) is saturated in the UV, u < 1. Substituting the 
wavefunction in the form (53) and using (56), (58) one obtains, 

5u 2 m) = k 2 ■ 2(i/ + l)B A 2+2u , (62a) 
r 2 ,2 1 x ' 2v B A 2 m 2u . . 

H w =* v^i))WA)r n£1, (62b) 

where 

B= (i^)'f^ 9 (f-l±f) { ^. (63) 

Note that we have exploited the fast falloff of the integrand at u > 1 and extended the 
integration to infinity. One observes that for m n ~ A the corrections to the velocities of 
the modes are suppressed by (A/ A*) 2+2u . For larger masses the corrections increase and at 
m n > A we obtain 

5u 2 nXl) ^k 2 BAm-^ . 



= k 2 ^ - tfj " / dxx^(J v {x)f , n > 1 . (60b) 

\Jv k ,/A Jo 



28 As mentioned in Sec. 3.1, this case is actually always realized in the model at hand, unless d>7. 
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Still this is parametrically small for all m n <C 1. 

The structure of the LV corrections to the dispersion relations found above fits naturally 
into the holographic interpretation of the model. The expressions (60), (62) are interpreted 
respectively as the contributions of the first and second operator in the Lagrangian (48) 
describing the perturbation from the LI fixed point. 



3.3.3 Higher-order corrections 

In principle, the quantum-mechanical perturbation theory allows to compute higher-order 
corrections to the dispersion relations of the KK modes. The Ith correction will involve 
sums over products of / matrix elements of the LV potential V(u), see (50). As the latter is 
proportional to k 2 , one will obtain in this way the dispersion relation in the form (19), as a 
series in powers of the spatial momentum. This is precisely the form used in the literature to 
analyze the phenomenological constraints on LV. However, the complexity of the formulas 
increases rapidly as one goes to higher orders /, making calculation of the coefficients in 
this series impractical beyond the case 1 = 2, which we are now going to consider. Not to 
overload the paper, we restrict the analysis to the zero mode. Apart from giving the scale M 4 
that suppresses the k A term in the dispersion relation, this analysis will provide information 
about the domain of validity of the perturbative calculation. 

The starting formula for the second-order correction to the frequency of the zero mode is 

oo 

K(2) = -E m n 2 |^l y l^)! 2 - ( 64 ) 

71=1 

Note that this is negative implying that the k A contribution comes with the minus sign. This 
appears to be a general feature of the holographic setup valid beyond the specific model of 
this paper. It is a direct consequence of the fact that the second-order correction to the 
frequency of the lowest KK mode is always negative. 

The evaluation of the matrix element appearing in (64) proceeds differently depending 
on which of the quantities, a or 1 + u, is larger. We have already encountered this situation 
above in the computation of the first-order LV contribution. Thus we have: 

Case (i): a < 1 + v. The integral defining the matrix element is saturated at « > 1, 
which yields, 

<0o|V10„> = fc 2 (2/ 00 - ™?, / dx x 1+V ~ 2a Ux) • (65) 

\Jv\ m nl -<V| JO 

We need the asymptotics of this expression at n ^> 1. There are again two possibilities 
depending on whether or not the remaining integral converges at m n /A — > oo. Using the 
formula 

t / x /2 / TXV 7I"\ 

J v [x) ~ \ — cos I x at x — >• oo 

y ' V irx \ 2 4/ 
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and integrating by parts we obtain the leading behavior, 

»m n /A 

„l+v-2a 



dx x l+u - Ia J u (x) = consti + const 2 • (m n / 'A)-**"-** . 

In this derivation we have used that the masses are given by the zeros of J u+ i, hence at 
n > 1, 

( 1XV 7T \ 

I IT + 4 + nn J 

In this way we obtain the following expression for the matrix element at large n, 

(0o|V|0„)ocA; 2 (2/ oo -Ol)A 2a x 



n~ 2 , a<-l/4 + i//2 

n-*-" +2a , a > -l/A + v/2 . 



Substituting this into (64) we observe that the sum over n converges provided a < l + v/2. 
To understand the origin of the divergence that appears in the opposite case, one recalls that 
Eq. (65) is valid only for masses smaller than the LV scale, m n <C 1. The divergence is an 
artifact of this low-mass approximation: the total correction to the energy level is of course 
finite in quantum mechanics. For a crude estimate we can just cut off the sum at n max ~ A -1 . 
Combining everything together we arrive at, 

K( 2 )--^-^) 2 x| A2+2 ^ a>1 + u/2 . (66) 

From this expression one infers the suppression scale M 4 . Restoring the explicit dependence 
on A* one obtains, 

M 4 ~ ( A - ( A " /A |;r ' a<l + U ' 2 (67) 
\A, (A./A) , a> l + v/2. K 1 

Note that M4 is parametrically higher than A* for a > 1/2. 

Case (ii): a > 1 + v. The matrix element is dominated by the region u < 1. Thus, we 
use Eqs. (53), (56), (58) and obtain, 

r(l + i/) \J u {m n /A)\ 

where B is defined in (63). Substitution of this expression into Eq. (64) produces a divergent 
sum, which we cut off at n max . The final result reads, 

5ul (2) oc -k A B 2 A 2+2v . (68) 

Note that the parametric dependence on the confinement scale A is the same as in the lower 
case of Eq. (66). Correspondingly, the mass parameter M 4 is given by the lower case of (67). 
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Let us discuss the conditions for the validity of the quantum-mechanical perturbation 
theory used above. A necessary requirement is that the correction to the eigenvalue is smaller 
than the spacing between the adjacent levels of the unperturbed Hamiltonian. Besides, the 
second-order correction must be smaller than the first one. This amounts to the conditions, 

< A ' > Suj lm < • ( 69 ) 

Clearly, they restrict the values of the spatial momentum k where the perturbative formulas 
make sense. By inspection of (60a), (62a) we find that the first inequality in (69) implies 

k < min{A 1_a , A~ u } , 

while Eqs. (66), (68) and the second inequality lead to 

{A 1 "" , a < 1 + v/2 
A- 1 -^" , 1 + v /2 < a < 1 + v 
1 a > 1 + v . 

One observes that the bound following from the second order of the perturbation theory is 
stronger for a > 1 + v/2. We conjecture that considering even higher orders will lower the 
bound further in the range 1 < a < 1 + z/, so that eventually one will arrive at 

k < minjA 1 " , 1}. (70) 

In other words, the upper bound on the spatial momentum for which the perturbative 
expansion of the dispersion relations works, is smaller or equal to the LV scale A*. This is 
quite similar to the situation in the RS-type model of Sec. 2. 



3.3.4 Zero-mode dispersion at large momenta 

One would like to go beyond the limitation (70) and study the dispersion relations at large 
momenta. In particular, it is interesting to understand if deviations from LI remain small 
when the momentum exceeds the LV scale A*. It turns out that, as far as the zero mode is 
concerned, this question can be easily answered in the affirmative if the parameter £ is zero 
or negative. Note that according to Eqs. (60) the first corrections to the velocities of the 
bound sates in this case are positive meaning that these modes are 'superluminal'. 

The argument is based on the variational theorem that states that the average of a 
Hamiltonian over any function is larger than the energy of the ground state. Consider first 
the total Hamiltonian + V. Its lowest eigenvalue gives the frequency of the zero 
mode. Thus we have, 

u 2 Q -k 2 < (<j>\H m + V\4>) 

for any function (f)(u). Let us substitute here the ground state eigenfunction O of the 
unperturbed Hamiltonian . Recalling that the lowest eigenvalue of the latter is zero we 
obtain, 

co 2 - k 2 < (<h\V\<h) = <H 2 (i) , (71) 
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where the r.h.s. has been computed in Eqs. (60a), (62a) for the cases a < 1 + v and a > 1 + u 
respectively. We conclude that the deviation of the (phase) velocity of the mode from unity 
is bounded from above by a small quantity behaving as a positive power of the IR scale A. 

To obtain a lower bound we use the property of the domain wall solution mentioned in 
Sec. 3.1 that f(u), j(u) are monotonically decreasing functions. For £ < this implies that 
the potential V(u) is also monotonically decreasing. Thus, it always exceeds 

Knin = fc 2 (2/ 00 -^)A 2Q , 

which is attained at u — A -1 . We now have the following chain of relations, 

"1 - k " = (CI#(o) + VW> > (<f>o X \H(o)W) + V mm > V mm , (72) 

where 4>q x is the exact ground state wavefunction of the full Hamiltonian + V. In the 
last inequality we have again applied the variational theorem, this time to the unperturbed 
Hamiltonian H(o). Eq. (72) gives the desired lower bound on the deviation of the mode 
velocity from one. We see that in the case under study (£ < 0) it is positive and proportional 
to A 2a . Combining the bounds (71) and (72) one concludes that for a large hierarchy between 
the IR scale A and the LV scale A* the dispersion relation of the lightest bound state is 
confined within a narrow wedge close to the lightcone. 

Interestingly, we see that the Lifshitz model has the same property as the Randall- 
Sundrum model of Sec. 2: the dispersion relation of the zero mode displays the emergent 
relativistic form at all momenta, including k ^> A*. Again, we observe that this property 
correlates with the superluminality of the bounds states. Similarly to the Randall-Sundrum 
case, this behaviour can be understood as due to the additional contribution to the potential 
(50) representing an effective k— dependent mass term peaked in the Lifshitz part of the 
geometry. In the superluminal case the squared mass is positive so that the zero mode 
is more and more repelled from the UV at large momenta, leading to the recovery of the 
emergent speed. The same reasoning suggests that the persistence of LI at all momenta must 
also hold for the excited KK modes, though we did not find a simple way to demonstrate 
this explicitly. To sum up, it seems rather generic that the emergence of LI is more efficient 
for superluminal models, at least within the validity of the holographic correspondence. 

4 Summary and Discussion 

In this paper we have studied the emergence of Lorentz invariance (LI) at low energies 
in strongly coupled Lorentz violating (LV) scale-invariant field theories (SFT's). We have 
considered two scenarios: relativistic CFT's which are perturbed at high energies by coupling 
to an 'external' LV sector; and Lifshitz SFT's which flow in the IR to a LI fixed point. 
Using explicit models that provide the holographic description of these scenarios, we have 
analyzed the implications of the emergent LI for observable quantities, such as correlators and 
dispersion relations of the bound states, and developed a systematic approach to estimate the 
LV corrections. We have confirmed that the strong dynamics accelerates the renormalization 
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group (RG) flow towards LI at low energies. In agreement with the general RG expectations, 
we found that the leading LV corrections are power-law suppressed by the ratio of the IR 
scale to the scale of LV, with the exponent related to the dimension of the least irrelevant 
Lorentz violating operator (LILVO). We have identified the IR scale that must be substituted 
in the RG formulas. For the correlation function in the Euclidean domain this is nothing 
but the inverse distance/time, at which the correlation is measured. For the bound state 
spectrum the relevant IR scale coincides with the scale of confinement; this can be much 
lower than the absolute energy or momentum of the particle. In other words, the leading 
corrections to the dispersion relations have the form of the shift in the propagation velocities, 
u 2 ~ (1 + 25c) k 2 , with 

5c oc (A/A*) A ~ d , 

where A* is the scale of high-energy LV, A is the confinement scale and A is the dimension of 
the LILVO. These findings are common to both scenarios that we considered and we believe 
them to be valid for any system where LI emerges due to strong dynamics. 

On the other hand, the nature and dimension of LILVO is model-dependent. In the 
Randall-Sundrum-type model realizing the first scenario the role of LILVO is played by 
the kinetic term <p 2 of the 'external' LV scalar field which couples to a CFT through a 
scalar operator 0$. Close to the IR fixed point acquires an anomalous dimension, so that 
dim0 = d — dimO^,. Then for the dimension of LILVO we obtain 

Ai = 2 + 2d -2 dim 0^ . (73) 

The case dimO^ — d/2, which is the minimal dimension of the scalar operator achievable 
in this setup, gives the most efficient recovery of LI with the leading corrections, such as 5c, 
suppressed by two powers of A/ A*. 

It might seem that a model with dimension of 0$ lower than d/2 (but still larger than 
the unitarity bound d/2 — 1) would lead to a larger A and hence stronger suppression of 
LV corrections. However, this is not true. Such a model will always contain a LV operator 
constructed of 0<f, itself, namely O 2 , whose dimension is 

A 2 = 2 + 2dimO . (74) 

This is smaller than (73) for dimO^ < d/2. Thus O 2 becomes LILVO and the suppression is 
actually weakened. This is precisely what happens for the alternative interpretation of the 
model of Sec. 2 mentioned in footnote 7. 

More constraints arise for a CFT that emerges as the IR endpoint of a Lifshitz flow. It 
contains the additional LV perturbation by the irrelevant spin-one operator, 

^ = ^ > (75) 

with a = dim(O^) — d > 0. The latter imprints a modified scaling in the observables. In 
particular, the LV correction in the speed of the scalar bound states receives an additional 
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contribution scaling like 29 5c ~ (A/A*) a , which competes with the intrinsic correction that 
scales according to (74). Whenever dimO^ is bigger than d/2 — 1 + a, the leading correction 
comes with the exponent 2a and becomes independent of dimO^,. In the Lifshitz flows of 
[26, 28] a cannot be too large 30 (a < 0.35 in d = 4) and therefore the LV effects are not very 
strongly suppressed. In view of this, it would be interesting to search for other (holographic) 
models with UV Lifshitz scaling where the dimension of LILVO would be higher. 

We have also analyzed the next-to-leading LV corrections to the dispersion relations - 

the contributions into uo 2 of the form k 2l M^~ 21 with / = 2,3, We have found that the 

mass scale M 2 i that suppresses the order-/ correction is given by A* multiplied by a certain 
power of the ratio A/A*. This power is related to the dimension of the LILVO in a non-trivial 
way; furthermore, it is different for the terms of different order /. In general, it can be both 
positive or negative depending on the model parameters, and correspondingly the scale M 2 ; 
can be lower of higher than A*. It is unclear at the moment if the scaling behavior of M 2 i 
can be determined from the general RG arguments or if it is strongly model-dependent. 

We have observed that the Taylor expansion of the dispersion relations in the powers of 
momentum inevitably breaks down at momenta at most as high as A*. This happens even if 
the individual terms in the expansion are suppressed. In the Randall-Sundrum type model of 
Sec. 2 it were possible to analyze the dispersion relations rather explicitly past this critical 
momentum and we found that, quite surprisingly, they tend to relativistic form at high 
k ^> A* if the 'external' LV sector is superluminal. In the subluminal case the recovery of 
the relativistic dispersion relations still holds for the excited bound states, while the lightest 
(massless) state gets non-relativistic at high momenta. For the model of the Lifshitz flow 
studied in Sec. 3 the dispersion relations cannot be determined analytically. Still, for the 
parameters in a certain range corresponding again to superluminal type of LV, we were able 
to obtain constraints on the frequency of the massless state which show that its dispersion 
relation lies in a narrow wedge close to the light-cone. In other words, its phase velocity u/k 
is constrained to be close to 1 at arbitrary k. It would be interesting to understand whether 
this property — persistence of almost LI dispersion relations at all momenta — is specific 
to the class of holographic models studied in this paper or is more generic. 

Our results have implications in two areas. First, the holographic models of the type 
studied in this paper may be relevant for the description of emergent LI in condensed matter 
systems. The phenomenon of emergent LI has been known for a while in the condensed 
matter literature, see e.g. [9, 10, 11], where it has been analyzed with a variety of methods 
including 1/N-, e-expansions and the exact RG. However, as we discussed, the emergence 
of LI becomes truly efficient only when the theory gets strongly coupled, which renders the 
traditional perturbative methods unreliable. The holographic duality provides a complemen- 
tary approach to assess the emergence of LI in strongly coupled condensed matter systems, 
which allows a systematic calculation of the LV corrections and their scaling exponents. Ad- 

29 The fact that the exponent equals 2a and not just a is a consequence of the additional discrete symmetry 
0^ A i \ —O^ appearing at the IR fixed point. Due to this symmetry the perturbation (75) contributes to the 
dispersion relations of the scalar bound states only at the second order of conformal perturbation theory. 

30 Othcrwise, the same bulk spin 1 field cannot give a relevant deformation of the UV Lifshitz point. 
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mittedly, an important challenge along these lines is to construct holographic models that 
would faithfully represent the essential properties of real materials. 

A related phenomenon, which also seems to be ubiquitous in condensed matter context, 
is emergence of rotational invariance (or isotropy). Despite the fact that any material lacks 
rotational invariance at the atomic/molecular level, many of them display an isotropic re- 
sponse in the continuum limit. Again, the approach to the isotropic fixed point can be 
described using the holographic methods by an analytic continuation of our discussion 31 . It 
would be interesting to establish the connection between the results of this paper and the 
previous studies of the emergence of isotropy in O(N) models [39] and lattice field theories 
[40]. In particular, it was found in these references that in the continuum limit all anisotropic 
corrections scale with the value of the exponent close to 2. It would be interesting to see 
what is the relation between this and our findings regarding the optimal LILVO dimension 
in the case of the Randall-Sundrum type model. 

Another important application of the mechanism for emergent LI studied in this paper 
appear in non-relativistic quantum gravity scenarios [2, 3]. As mentioned in the Introduction, 
the idea (as in technicolor/compositeness models) is to embed the Standard Model (SM) 
fields as the lightest excitations of some strongly coupled sector which has order-one LV 
in the UV. The strong dynamics would then enforce the approximate LI seen in the known 
particle physics. This requires realizing the SM fields as (partially- or completely-) composite 
states. The most minimal holographic realization of a completely composite SM in this 
context is the Randall-Sundrum I model [25], with the entire SM living on the IR brane and 
LV physics restricted to the UV brane. In this model there is no tower of states with the 
quantum numbers of the SM fields, so at leading order in 1/N there is no LV whatsoever 
in the SM. Only the inescapable tower of spin 2 resonances receive small LV corrections at 
low momenta in this approximation. Of course, quantum effects in the bulk are expected 
introduce some 1/iV-suppressed LV into the SM sector. We leave the interesting question of 
estimating these effects for future. 

In order to be more generic, one can place the SM fields in the bulk, realizing in this 
way 'partial compositeness'. This will require to generalize the construction of the present 
paper to include fermions and gauge fields. While inclusion of fermions is expected to be 
straightforward and closely follow the track presented in this paper for scalars 32 , one expects 
some qualitative differences to arise for the gauge fields. From the bulk perspective, the 
reason is that the gauge-field zero-modes are delocalized in the bulk: their wave functions 
penetrate into the UV and so they will be sensitive to the LV UV physics. In the CFT 
language the gauge field in the bulk is dual to an 'elementary' 4-dimensional gauge field 33 
coupled to the CFT through a conserved current with protected dimesion, dim = 3. Such 

31 We actually performed such continuation when we studied the correlation functions. 

32 Fermions can be given different degree of compositeness by adjusting the bulk and brane mass terms [41]. 
In the present context, since the deviations from LI are most tightly constrained for the first generation of 
quarks and leptons, these are the fermions that must be the most composite. 

33 To give rise to a zero-mode the bulk gauge field must obey Neumann boundary conditions on the UV 
brane. Its boundary value is then interpreted as the 'elementary' field [42]. 
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coupling is marginal and implies logarithmic sensitivity to the UV. We have observed this 
kind of behavior in the model of Sec. 2 when the 'elementary' scalar couples to the CFT via 
operator of dimension d — 1 (see the end of Sees. 2.1, 2.2). In this case the LV corrections 
are indeed suppressed only logarithmically. This represents a serious obstruction, since it 
leads to the reappearance of the Lorentz fine-tuning problem. 

Let us point out a possible way to circumvent this problem. It relies on the observation 
that even though the LV corrections log-persist along the RG flow, they are inversely pro- 
portional to the coupling constant. Namely, the logarithms enter into the formulas in the 
combination g^ 2 (log A^/A) -1 with g 4 ~ \/L/b being the coupling between the 4-dimensional 
gauge field and the CFT (cf. Eqs. (11), (20)). This suggests that there should still be a 
limit where the LV corrections are suppressed, involving again strong coupling, g± 1. In 
hindsight, one realizes that this limit corresponds to enhancing the gauge kinetic term in the 
IR part of the geometry. A dynamical realization of this might be possible by introducing 
interaction of the bulk gauge field with a dilaton, see [43]. In this case the 4-dimensional 
gauge boson would become composite and the gauge symmetry would be itself emergent. 
This scenario deserves further study and is left for future investigation. 

Similarly, we can foresee what happens to gravity in this framework by looking at the 
model of Sec. 2 in the case when the scalar operator has dimension 4 (for d = 4). Since 
this leads to an irrelevant coupling, LI does not emerge in this sector. This property may 
be actually an advantage. It implies that it is possible to have sizable LV in gravity without 
destroying the stability of LI in the matter sector. This is interesting from the phenomeno- 
logical viewpoint as it opens an opportunity to probe gravitational LV. The existing data 
already place non-trivial bounds on LV in gravity [47]. However, these bounds are much 
milder than in the matter sector, and so only a moderate amount of tuning is required to 
satisfy them. 

Finally, let us emphasize that the existence of the strong-coupling mechanism for emer- 
gence of LI together with the previously known ones (based on a non-relativistic SUSY 
[22] and on a separation of scales [44]) opens up an interesting connection between parti- 
cle physics phenomenology and non-relativistic gravity. Generally, avoiding the fine-tuning 
problems associated with fundamental LV will require some mechanism — some 'new physics' 
- operating down to a very low energy scale. Then via naturalness the non-relativistic grav- 
ity models become very predictive. A very rough estimate of the typical new-physics scale is 
around 10 orders of magnitude below the LV scale. The latter is bounded from above by the 
gravitational observations. Taking it to be ~ 10 15 GeV [45, 46], one obtains a new physics 
scale as low as ~ 100 TeV. By improving the constraints on the new physics scale on one 
hand and on the LV scale on the other, one may eventually cover the whole window required 
for the mechanisms of emergent LI to operate. In our opinion, this further motivates the 
experimental searches for new physics, including Lorentz Violation. 
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A Bound states in RS model with irrelevant coupling 

Here we consider the model of Sec. 2 with the choice of parameters such that v > 1. Ac- 
cording to the discussion in the main text, in this case one does not expect appearance of 
LI at low energies. We confirm this expectation explicitly by studying the spectrum of KK 
modes. 

The KK spectrum is still determined by the basic Eq. (14). Let us restrict to the choice 
Huv — — and make a technical assumption that v is non-integer. This will simplify ma- 
nipulations with the Bessel functions. Expanding the contributions with the "UV" -subscript 
at pL <C 1 one obtains from (14), 

(1 + x)u 2 - (c 2 + x)k 2 ( v\ 2v 1 , /A x 4r(2-z/ 



xA 2 



(^) " ViCp/A) = ^^{LKf^J^ip/K) , (76) 



where x has been defined in (9). The r.h.s. here is suppressed by the small factor (LA) 2( - U 1 \ 
Neglecting this contribution one obtains a LV gapless mode with the dispersion relation 

^ = c l±^l k 2 (77) 
1 + x 

and a family of LI massive modes, 

«J = k 2 + (AjS) 2 , n = l,2,..., (78) 

where are positive roots of J v -\. The formula (77) clearly matches with the two-point 
function of 0, Eq. (8), and thus the corresponding mode is identified with the elementary LV 
scalar whose velocity has been renormalized by the interaction with the CFT. On the other 
hand, the family (78) gives just the standard KK spectrum of the bulk scalar in AdS with 
the Dirichlet boundary conditions on the UV brane. So from the dual viewpoint these are 
the bound states of the CFT decoupled from the LV deformation. 

It remains to see how this picture is affected by taking into account the r.h.s. of (76). 
One finds that the frequencies of the CFT modes get shifted by 

^ 2=A2 2F(2 - y) (jSr+V^S) x(LA) 2 ^ 

r W Jl-i(&) (l-c 2 )A; 2 + (l + x)(A J t ) 1 ) 2 ' 
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Figure 4: Dispersion relations of the eigenmodes in the leading approximation (a) and after 
taking into account the mixing (b). 

For c 2 < 1 this is basically the end of the story: the corrections (79) are uniformly small 
at all spatial momenta, so the spectrum indeed splits into almost LI CFT plus a single LV 
mode. 

However, if c 2 > 1 the corrections blow up at 



implying that close to these points the perturbative scheme breaks down. To understand 
what happens one observes that (80) corresponds to the intersections between the dispersion 
relations (77) and (78). Thus we have a situation analogous to the level- crossing in quantum 
mechanics: as far as we completely neglect the r.h.s. in (76) there are two eigenvalues (fre- 
quencies) that cross at a given value of spatial momentum. It is known that the perturbation 
actually prevents the levels from crossing, so that the dispersion relations of the modes are 
reconnected, see Fig. 4. This produces a rather peculiar situation when the LV mode (77) is 
"cut in pieces" and its chunks are inserted between the LI KK modes. As a result we obtain 
a tower of modes, all of which are essentially LI at large momentum, while the massive modes 
are also LI at small k. However, every mode strongly violates Lorentz invariance in some 
range of momenta (for the gapless mode this range includes k = 0), so this case clearly does 
not exhibit emergent LI. 

B Lifshitz flow: Corrections to the correlator in IR 

Un this Appendix we emplement the perturbative scheme described in Sec. 3.2 to determine 
the LV corrections to the correlator (46). The reader should refer to that section for the 
notations and conventions. 




(a£\) 2 



(80) 
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B.l Solutions in the a'-region 

We start with the a'-region where we need to find the solutions to Eq. (37) up to the order 
0(w 2 ,k 2 ). Treating the terms containing w 2 ,k 2 in (37) as a perturbation, one uses the 
standard "variation of constants" technique and writes, 

C » = P aT (u)^l(u) , (81) 

where the indices a, r take the values ± and the coefficient functions are required to satisfy, 



This yields, 



P-{u) = 1- T fmtt^X , P- + (u) = [ U fnu){<PtY , (82a) 

P + _(«) = -£ d ^K,{u) (0i°l) 2 , P ++ (u) = 1 + £ ^/C(^°i0i°l , (82b) 



IC(u)=u^g[ 1 -±^w 2 + fk 2 ) , 



where 

K.(n\ = i t l - d n(— j 

and Mo is an arbitrary normalization point in the c-region, i.e. 1 <C u <C w~ l ,k~ l . In 
deriving (82) we have used the expression for the Jacobian of the functions 0^,°-? 

that follows from their asymptotics (43). 

Similarly, for the solutions defined by the asymptotics (38) we obtain, 

k-(u) = *£><«) (i - £ f Jf-\_) K ^) + «(») £ f Jf-\_f «- ) 2 • 

(83a) 
(83b) 

where we used the Jacobian 

do) AO) i _ /o(A+ - A_) gu d - x 



9- f 



Note that all but one integrals appearing above converge at u — > which allowed to extend 
the domain of integration down to this point. On the other hand, the integral entering the 
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second term in (83a) diverges at u — > 0, so we have normalized it at an arbitrarily chosen 
point «] < 1. We will see below that the final answer does not depend on this choice. 
Combining (83) with (81), (82) one arrives at the analog of the relations (45), 



c , CT = (v; { T o) +x CTT )0 a , 7 



with 



f u ° g.du 



«.(°) 

) 



_ f uo g^du (Q) (o) _ (o) [ U1 g-du (Q) (o) 

r° 9-du r (o) , (0) 

_ r° p-rfix (o) (o) (o) K AQ) AQ) 

X++ ~ " L /o(A + " X-f^'- V++ Jo JojK^Xl) ' ,+ ' 

In deriving these expressions we made use of the formula for the determinant of the ma- 



trix Vi% 



\V®\ = 



W[$U®] /o(A + -A_)" 



B.2 Solution in the 6- region 

In the 6-region one has to take into account the deviation of the background functions from 
the constant values, see Eqs. (34). The solution is again found by the "variation of constants" : 

<f )b (u) = R(u)^\u) + Q(u)^l(u) , 
where (f>^ is given by (41) and 

is the second linearly independent solution of the unperturbed equation. The coefficient 
functions must satisfy 

R'^T + Q'^i = o , R'<t>f + Q'^l' = Hu) , 

where 
This yields 

R(u) = l- duu^Fiu^liu) , Q(u)= duu 1 - d T(u)<l> { °\u) (84) 

J UQ J OO 
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where we have used the Jacobian 

Note that the integral in the expression for Q(u) is taken from infinity to ensure that the 
solution vanishes at u — > oo. 



B.3 Matching solutions in the c-region 

We need to relate the solution (fib to (fi c ,±, 

cj) b (u) = U-(j) c -(u) + U + (j) c>+ (u) , 

where U± are constants. To this end we expand (fib, Cj ± in the c-region in the powers of u 
and match the corresponding coefficients. It will be enough to keep track only of the terms 
containing u d ^ 2±u . 

We start with the expansion of 4> c> ±. To evaluate the integrals entering (82) one uses the 
asymptotic expressions at u >> 1, 



K{u) ~ ^"'[(l + ( 9oo + Uu- 2a )p 2 E + (til - 2f 00 )u~ 2a w 2 } , 
u d/2±"(i + a± u~ 2a ) , a± 



c,± I 1 + a ± M J > a ± - A , ^ x • 

Note that one should be careful to take into account the subleading terms arising from the 
tails of the domain wall background. Substitution into (82), (81) yields, 

<fc,- = « d/2 -"{i + ^ + 4l/ "°T a) + /- + «+ + + (Ot> - 2/«>- 2 ] } 

//2+ 1 if 5 ! + sf^)^ + + 2a -^ + (f - - ^ } + • • • • 

(85a) 

+ i^r^ [<*» + /- + 2 ^ + - 2 ^ 2 ] } 

+ « d/2+ "{l - ^ - i ^f-^y foco + /oo + « + + a_)p| + (£& - 2/o> 2 ] } + ..., 

(85b) 

where dots stand for the terms with other powers of u. 

To find the relevant terms in the expansion of the function (fib one makes two observations. 
First, the constant shift of the function R(u), see (84), amounts to the change of the overall 
normalization of (fib and thus will not affect the result for the correlator (which is of course 
independent of this normalization). It is straightforward to see that only the constant part 



u 



u d/2-u < 



41 



of R(u) contributes into the coefficients of interest. We will conveniently normalize this 
contribution to f . Second, the function Q(u) can be evaluated explicitly using the formula, 



J V 



j -i-2a( T ^ f n\2 ^T(-a)Y(-a - v)T(-a + v) 
dxx 1 2a (K u (x)) = 



4r(l/2-a) 



ny 2a . i 

2 F 3 -, -a; 1 - a, 1 - v, 1 + v; y 



4z/asin7rz/ \2 

cg,v -2a.-2v / j 

+ -rr r- (I») 2 2 F 3 - - z/, -a - z/; 1 - a - z/, 1 - 2z/, 1 - z/; r 



(a + z/) 

2~2u -2a+2v / j \ 

+ 8 (a-i/) (r( ~^ )2 2Fs ( 2 + + ^ 1 ~ " + l/ ' 1 + 2Z/ ' 1 + Z/ ' y J ' 

where 2-F3 is the generalized hypergeometric function. It is an entire function of its last 
argument and thus adds only integer powers of y 2 . Given this information, we obtain, 



0, = IV/ 2 -" + f/jV/ 2 +' [l + s a (2 9ooP 2 E + (til - 2/ oc » s 2 + 2 « 

+ s 1+a (a(d + 2a)(f OQ - g^) + 2g O0 (nL) 2 )p 2 E Jl 

where are given in (44) and 

sin7rz/r(l — a)r(l — a — u)Y(\ — a + v) 



(86) 



Sty 



2^(3/2 - a) 
Combining (86), (85) we find, 

( „.2J2 r,,2-2oL 

U- =U® 1 1 - ^ - ^--^ [{900 + /oo + «+ + a_)p| + - 2 /oo ) W 2 ] 



2+2i/9 „,2+2i/-2a 



+ s a (2 gooP 2 E + (til - 2f OQ )w 2 )p E 2+2a + si +a {a(d + 2a)(f OQ - 9oo ) + 2 goo (fiL) 2 )p 2 E a 
B.4 Evaluation of the correlator 

We are now ready to evaluate the coefficients T± entering (39). A simple algebra yields, 

T a = Ui \V T f+l Ta ), 
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where keeping only the LV contributions, 



X_ 



w 



X, = w z 
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2-2a 



4i/(l -a) 



-V 



rf M /c^ |^ + (6l-2/ oo ) 
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,2-2a 



+- 



4z/(i - «; 



(0) "0 



2-2u-2a 



4z/(l - 1/ - a) 



- V_ 



(0) «0 



"4i/(1 + z/ - a) 
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(87a) 
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(87b) 



X-+ = ^ 2 - 



"1 



^0 



- - 2/00) 



n, 2 ~ 2a 

y(0) U 



n 2-2v-2a 
T/ (0) — ^ 

Av(l-a) ++ Av(l-v-a) 



+ 



w < — 



'c,+~Ta,- 



++ 
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++ 4z/(l -a) 



2+2z/-2a 

-+ 4z/(l + v-a) 



2+2a 



(87c) 
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(87d) 



and 



K{u) = 



g-u 



l-d 



1 + a 2 



f • 



/o(A+ - A_) " V / 

Note that expressing 01°^- via <p^l by inverting the relations (45) one can check that X CTT 
are independent of the normalization point uq within the accuracy of our approximation, as 
it should be. On the other hand, the dependence on u±, that affects the definition of the 
solution (p a -, is still present; it will drop out only at the last step of the calculation. 

Finally, we must compute the ratio T+/T1 that according to Eq. (40) will give us the 
correlator of interest. Expanding this ratio to quadratic order in /U^ and to linear 
order in X aT and dropping the local contributions polynomial in w 2 and k 2 we find, 



Qs oc 



c/i°) (vl°l) 2 



1 + 



I/(o) 



v^x ++ - k^x + _ - vl°li- + - ^x__ + 



2Vl°M°l 



(o) 



x_ 



u 



(o) 



jj(o) |y(o)| 



2\/ (0 V (0) 
z v ++ v +-.j 



(0) 



+ 



2vi°M°2 



(0) 



-+- 



(o) 



(Vi u J) 



(0)\2 



43 



Inserting here Eqs. (87) after a somewhat long but straightforward calculation, where one 
has to make use of the linear relations between the functions <f)^i and <f>^l, one obtains the 
expression (47) from the main text with 

sin7rz/r(l — a)T(l — a — u)T(l — a + v) 



Ci =(UL - 2/co) 

2 



20Fr(3/2 - a) 



(88a) 



V^lJo /o(A+-A_)" 



C,=2 



+ (Ul - 2/c 



2-2a 




/ 



_2i/(l - a) v[°l M 1 -v-a) 



(88b) 



-2i/- 
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(0) 
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Uo 



1-d 



(yi°2) 2 io /o(A + -A_)"V / 



1+0' 2 



(0) X2 



2+2^-2a 




+ 



4i/(l + i/-a) vd^M 1- ") VV^/ 4i/(l-i/-a) 



(0)\ 2 



If 



2-2^-2a 




where we have restored the explicit dependence on the functions defining the domain wall 
background. Note that these expressions are clearly independent of the arbitrary choice 
of the normalization point u\. As already noted before, despite the fact that Eqs. (88b), 
(88c) apparently contain the other normalization point Mo, they are in fact independent of it 
within the accuracy adopted in our calculation as long qsuq^I: the boundary terms from 
the integrals are canclelled by the power-law counterterms. Moreover, the corresponding 
contributions in the correlator (47) are actually important only when they dominate over 
the Ci-term. For the CVterm this happens if a > 1, while for the C^-term — if a > 1 + In 
these cases the integrals in (88b), (88c) become convergent and can be extended to infinity 
so that the counterterms disappear. 



C Energy- momentum of the IR brane 

In Sec. 3.3 we consider a brane that cuts off the large-w portion of the spacetime interpolating 
between the Lifshitz region and AdS. Here we show that the energy-momentum of such a 
brane can be decomposed into a constant negative tension plus a contribution satisfying 
the null energy condition (NEC). We will assume that the brane is located at an orbifold 
fixed point under a Z 2 reflection symmetry, similarly to the RS construction [25]. Then the 
energy-momentum of the brane is given by the junction condition (see e.g. [48]), 

where 7 M „ is the metric induced on the brane and K^ v is the brane's extrinsic curvature. 
The latter is defined as 
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where n u is the inward pointing (n u < 0) unit normal to the brane. A straightforward 
calculation using the form (22) of the metric yields 



Tr 



brane 
00 



d-l Lf 



AnG 



u 2 g 



Tl 



brane 



AnG 



(d-l)L Lf 
ugf 



u 2 g 



(89) 



where both expressions must be evaluated at the position of the brane u = A" 1 . We see 
that the energy density of the brane is negative. This is the common property of the IR 
branes in the context of the RS-type models. It can be accommodated by assuming that the 
brane has negative tension. Due to the orbifold symmetry accross the brane this does not 
lead to instabilities. In our case one can take any value of the brane tension satisfying the 
inequality, 

d-l 



a < 



AnGLg 



u=\- 



Subtracting the corresponding contribution from (89) one obtains the residual energy-momentum 
tensor T^ arie , such that for any ((d + l)-dimensional) null vector l M 



rpbrane ^> 



AnG ugf 



u=A- 



As discussed in Sec. 3.1, f{u) is a monotonically decreasing function on the domain wall 
solution implying that its derivative is negative. 
NEC. 



Thus we conclude that T^ v ane satisfies 
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